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We calculate and resum a perturbative expansion of nuclear enhanced power corrections to the 
structure functions measured in deeply inelastic scattering of leptons on a nuclear target. Our 
results for the Bjorken x-, Q 2 - and yl-dependence of nuclear shadowing in F 2 (x, Q 2 ) and the nuclear 
modifications to F^(x,Q 2 ), obtained in terms of the QCD factorization approach, are consistent 
with the existing data. We demonstrate that the low-Q 2 behavior of these data and the measured 
large longitudinal structure function point to a critical role for the power corrections when compared 
to other theoretical approaches. 
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In order to understand the overwhelming data from the 
Relativistic Heavy Ion Collider (RHIC) and make predic- 
tions for the future Electron Ion Collider (EIC) and Large 
Hadron Collider (LHC) in terms of the successful pertur- 
bative QCD factorization approach we need precise 
information of the nuclear parton distribution functions 
(nPDFs) 4>f{ x i M 2 ) of flavor /, atomic weight A, momen- 
tum fraction x and factorization scale \x. Although the 
/i-dependence of leading twist nPDFs can be calculated 
in terms of pQCD evolution equations, the x- and in- 
dependent boundary condition at a scale no must be 
fixed from existing measurements, mostly of the struc- 
ture function F^(x, /x 2 ) in lepton- nucleus deeply inelastic 
scattering (DIS) 2). It was pointed out recently Q that 
the available fixed-target data contain significant higher 
twist effects hindering the extraction of the nuclear par- 
ton distributions. 

On the other hand, it has been argued that for physi- 
cal processes where the effective x is very small and the 
typical momentum exchange of the collision Q ~ \i is not 
large the number of soft partons in a nucleus may sat- 
urate 0. Qualitatively, the unknown boundary of this 
novel regime in (x, Q) is where the conventional pertur- 
bative QCD factorization approach should fail 

In this Letter we present a pQCD calculation of the 
resummed power corrections to the DIS nuclear structure 
functions and demonstrate their importance in extracting 
nPDFs. From our results we quantitatively identify the 
characteristic scale of these power corrections £ 2 <C rn 2 N , 
with nucleon mass m^v = 0.94 GeV. We conclude that 
for Q 2 > m 2 N inclusive lepton-nucleus DIS is above the 
saturation boundary and can be treated systematically in 
the framework of the pQCD factorization approach with 
resummed high twist contributions. 

Under the approximation of one-photon exchange, 
the lepton- hadron DIS cross section dagh/dx dQ 2 cx 
L^W^(x,Q 2 ), with Bjorken variable x = Q 2 /{2p ■ q) 
and virtual photon's invariant mass q 2 — —Q 2 . The 
leptonic tensor L M „ and hadronic tensor W^ v are de- 
fined in 0. The hadronic tensor can be expressed in 
terms of structure functions based on the polarization 
states of the exchange virtual photon: W^ v (x,Q 2 ) = 



in 0- The transverse and longitudinal structure func- 
tions are related to the standard DIS structure functions 
Fi, F 2 as follows: F T (x,Q 2 ) = F^Q 2 ), F L (x,Q 2 ) = 
F 2 (x,Q 2 )/(2x) - Fx(x,Q 2 ) if 4x 2 m% < Q 2 g. In DIS 
the exchange photon 7* of virtuality Q 2 and energy 
v = Q 2 / (2xmpf) probes an effective volume of trans- 
verse area 1/Q 2 and longitudinal extent Azn x xjqjx, 
where Azn is the nucleon size, xn — l/(2romjv) ~ 0.1 
and ro ~ 1.2 fm |5(. When Bjorken x <C xjv the lepton- 
nucleus DIS covers several nucleons in longitudinal direc- 
tion while it is localized in the transverse plane. Although 
a hard interaction involving more than one nucleon is 
suppressed by powers of l/Q 2 @, it is amplified by the 
nuclear size if x <C xn- In the collinear factorization 
approach we evaluate the nuclear enhanced power cor- 
rections to the structure functions to any power in the 
quantity (£ 2 /Q 2 )(A 1 ^ 3 — 1) and keep £ 2 to the leading 
order in a s . 

In terms of collinear QCD factorization, the structure 
functions at the lowest order in a s are given by [|j 

4 LT) (x,Q 2 ) = ±Y,4tM*>&) + 0(a B ), (1) 
/ 

= 0(a s ) , (2) 



where (LT) indicates the leading twist contribution, 
runs over the (anti)quark flavors, Qf is their fractional 
charge and 0/ is the leading twist quark distribution: 



cj>f{x,ii 2 
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(p\Mo)^M*o)\p) (3) 



4" F T (x,Q 2 ) + ef F L (x,Q 2 ) 



where e^", 



fl IS 



Sj, are given 



in the lightcone A + = n^A^ = gauge for hadron mo- 
mentum p^ = p + h fl , where = [1,0,0 J and n M = 
[0,1, 0j_] specify the "+" and "— " lightcone directions, 
respectively. In Eq. 10 the parameter Ao = P + Uq ■ 

To compute the nuclear enhanced power corrections we 
choose the A + = gauge and a frame in which the nu- 
cleus is moving in the "+" direction, p M = P't/A = p + n fJ ', 
and the exchange virtual photon momentum is q^ = 
—xp + n fl + Q 2 1 (2xp + )n fJ - . In this frame the struck quark 
propagates along the "— " direction and interacts with 
the "remnants" of the nucleus. The leading tree level 
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FIG. 1: (a) Tree level multiple final state scattering of the 
struck quark in DIS on nuclei; (b) and (c) leading coupling 
for transversely and longitudinally polarized photons to the 
active partons, respectively; (d) effective scalar vertex for the 
basic two-giuon interactions. 



contributions to in lepton-nucleus DIS are given by 
the Feynman diagrams in Fig. ^a). The cut-line rep- 
resents the final state quark and the blob connected to 
7* - the lowest order coupling between the photon and 
the active partons. For transversely polarized photons 
Fig. ^b) gives the leading twist contribution in Eq. . 
For longitudinally polarized 7* the leading tree coupling 
Fig. Wic) results in the first power correction to F^: 

F L (x,Q 2 ) = F( LT '(*,Q 2 ) + i£Q2 4 ^0 
x / ^e^(p|^(0)^V/(A )F A 2 \p), (4) 



2p^ 



where the 0(a s ) F^ T ^ is given in |£| and the operator 



Ft = 



d~X n dX F+ a (\n)F+(\ ) 



2tt 



(P 



0(A„)fl(A o -A o ) ■ (5) 



Since the effective 7£ coupling is made of short-distance 
contact terms of quark propagators 9], the integrals of 
Ao and A„ in Eq. (JSJ should be localized. 

To derive the gauge invariant and nuclear enhanced 
power corrections at the tree level, we add gluon interac- 
tions to the struck quark and convert the gluon field op- 
erators in the hadronic matrix element of to the cor- 
responding field strength 0, QJ • We note that there must 
be an even number of interactions between the 7* cou- 
pling and the final-state cut because the quark propaga- 
tor of momentum Xip + q has only two types of terms |9j , 
i{"f + /2p + )/(xi— x±ie) (pole term -x->) and i(xp + /Q 2 )j~ 
(contact term —)—>), and the gluons are transversely po- 
larized. Consequently, the summation of the coherent 
multiple scattering can be achieved by sequential inser- 
tion of a basic unit consisting of a pair of gluon interac- 
tions, left-side of Fig.^d), connected by a quark contact 
term, and a quark pole term to the left (L) or right (R) 
if the unit is to the left or right of the final-state cut. 
Integrating over the loop momentum fraction Sj , we can 
replace this unit by an effective scalar interaction, right- 
side of Fig.^d), with a rule, 
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2tt Xi — Xi-i — it 
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.ziElhl R 



Xi—x—te 



(6) 

In Eq. |JBJ the boost-invariant operator F 2 (Xi) is different 
from in Eq. JSJ and is defined as 



dXi F+ a (X i )F+(X i ) 



2tt 



0(X> - Xt) 



(7) 



Its expectation value can be related to the small- 
x limit of the gluon distribution, (p\F 2 (Xi)\p) s» 
lirn^o \ x G(x, Q 2 ), and is independent of A^. It is the 
J dXi in Eq. 10 that gives the A 1 / 3 -type nuclear enhance- 
ment [f| . Assuming that the interaction leaves the nucle- 
ons in a color singlet state, taking all possible final state 
cuts [T?il | and performing the integrals over the incoming 
partons' momentum fractions we obtain the twist 2 + 2n 
contribution to the structure functions: 
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— (P A $ f (0)^^X0) F 2 J] yj dA^(A,)F 2 (A 



(8) 

\Pa) , (9) 



with corrections down by powers of the nuclear size. © in a model of a nucleus of constant lab frame density 

We evaluate the multi-field matrix elements in Eqs. ©, P( r ) = 3 /( 4 ™o) and approximate the expectation value 



of the product of operators to be a product of expectation 
values of the basic operator units in a nucleon state of 
momentum p = Pa/ A: 



n n 

(p A \ do \p A ) = 4 \p) n [ n p (p\Oi\p) 



with the normalization N p = 3/ (8irrQmp{). The integrals 
J d\i9(Xi) — (3romAr/4)( J 4 1 / 3 — 1) are taken such that 
the nuclear effect vanishes for A = 1. Resumming the 
A 1 / 3 -enhanced power corrections Eqs. JSJ, © we find: 



F*(x,Q 2 
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£2 (A l/3 _ X ) 



„d»4 LT) (x,Q 2 ) 



iff 1 I , 



F?(x,Q 2 ) a AF^ T \x, 



Q 2 



(10) 



^2(^1/3 _ X) 



iV 
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i /4£ 2 



n! V Q 2 



d"4 LT) (x,Q 2 ) 



AFf T) (x,Q 2 ) + i^F^(a;, Q 2 ) , (11) 



where iV is the upper limit on the number of quark- 
nucleon interactions and £ 2 represents the characteristic 
scale of quark-initiated power corrections to the leading 
order in a s 

e = ^0^( P \P 2 { x im . 

° r o 

In deriving Eqs. itTUI) . (fTT|> we have taken {p\ F 2 q \p) w 

(3romAr/4)(p| F 2 {\i) \p) and N ss 00 because the effec- 
tive value of £ 2 is relatively small, as shown below. 

Eqs. HI Uf l. (|ll|l are the main result of this Letter. Im- 
portant applications to other QCD processes and ob- 
servables that naturally follow from this new approach 
are given in The overall factor A takes into ac- 

count the leading dependence on the atomic weight and 
the isospin average over the protons and neutrons in the 
nucleus is implicit. We emphasize the simplicity of the 
end result, which amounts to a shift of the Bjorken x 
by Ax = x$ 2 (A 1 ^ 3 — 1)/Q 2 with only one parameter 
£ 2 oc lim^—^o xG(x, Q 2 ). In the following numerical eval- 
uation we use the lowest order CTEQ6 PDFs [T^ |. 

Fig. |21 shows a point by point in (x, Q 2 ) calculation 
of the process dependent modification to F2(A)/F2(D) 
(per nucleon) in the shadowing x < 0.1 region compared 
to NA37 and E665 data [HQ. We find that a value 
of £ 2 = 0.09 - 0.12 GeV 2 , which is compatible with the 
range from previous analysis |l5j of Drell-Yan transverse 
momentum broadening (£ 2 ~ 0.04 GeV 2 ) and momen- 
tum imbalance in dijet photoproduction (£ 2 ~ 0.2 GeV 2 ), 
makes our calculations consistent with the both x- and 
A-dependence of the data. Our calculations might have 
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FIG. 2: All-twist resummed F 2 (A) /A (Z3) calculation from 
Eqs. JT2J|, JTIJ versus CERN-NA37 O and FNAL-E665 
data on DIS on nuclei. The band corresponds to the choice 
£ 2 = 0.09 - 0.12 GeV 2 . Data- Theory, where Ad-t is com- 
puted for the set presented by circles, also shows comparison 
to the EKS98 scale-dependent shadowing parametrization 0. 



overestimated the shift in the region x close to xn where 
the TV ps 00 should fail [H. In Fig. we impose 
Q 2 = m 2 N for virtualities smaller than the nucleon mass, 
below which high order corrections in a s (Q) need to 
be included and the conventional factorization approach 
might not be valid. Our result is comparable to the 
EKS98 scale-dependent parametrization of existing 
data on the nuclear modification to F 2 A {x, Q 2 ), as seen in 
the Ad-t — Data — Theory panels of Fig. [3 We empha- 
size, however, that the physical interpretation is different: 
in [2( the effect is attributed to the modification of the 
input parton distributions at fio = 1.5 GeV in a nucleus 
and its subsequent leading twist scale dependence. In 
contrast, our resummed QCD power corrections to the 
structure functions systematically cover higher twist for 
all values of Q > fi . 

With £ 2 fixed, Fig. shows the predicted Q 2 depen- 
dence of F 2 (Sn) I F 2 (C) . The Q 2 behavior of our result, 
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FIG. 3: CERN-NA37 data EH on F 2 (Sn)/F 2 (C) show evi- 
dence for a power-law in 1/Q 2 behavior consistent with the 
all-twist resummed calculation. Comparison to Gribov-Regge 
model by Armesto et al. [TfH and the EKS98 shadowing 
parametrization 0] (in Ad-t) is also presented. The bottom 
right panel illustrates the role of higher twist contributions 
to Fl on the example of R = (tl/ct- Data is from SLAC- 
E143 0. 



distinctly different from a Gribov-Regge model calcu- 
lation by Armesto et al. [l^, highlights the important 
role of resummed QCD power corrections, when pre- 
sented versus NA37 data [T^ • Data-Theory panels again 



show comparison to the EKS98 parameterization Q . The 
QCD power corrections at small and moderate Q 2 to the 
leading twist structure functions can also be tested via 
the ratio of the cross sections for longitudinally and trans- 
versely polarized virtual photons 



R{x,Q 2 ;A) = 



<7j- 



F*(x,Q 2 ) 



(12) 



Comparison to E143 data for 12 C ll8| is given in the bot- 
tom right panel of Fig- EI T ne dashed curve includes only 
the bremsstrahlung and gluon splitting from Q 

and is insensitive to modifications of the nuclear parton 
distributions. The gray band represents a calculation of 
R from Eqs. ifTfl )) . lfTT ]l . 

In conclusion, for £ 2 = 0.09 — 0.12 GeV 2 our calculated 
and resummed power corrections are consistent with the 
x-, Q 2 - and A-dependence of existing data on the small- 
x nuclear structure functions without any leading twist 
shadowing. Our results, therefore, give an upper limit for 
the characteristic scale £ 2 of these power corrections in 
DIS on cold nuclear matter. Furthermore, the enhance- 
ment of the longitudinal structure function F^(x,Q 2 ) 
strongly favors a critical role for the higher twist in the 
presently accessible (small x < 0.1, Q 2 ~ few GeV 2 ) kine- 
matic regime. Less leading twist shadowing and corre- 
spondingly less antishadowing than currently anticipated 
will have an important impact on the interpretation of 
the d + A and Au + Au data from RHIC 19]. The 
predictions of this systematic approach to the process- 
dependent low Q 2 nuclear modification in QCD processes 
will also be soon confronted by copious new data from 
RHIC, EIC and the LHC. 
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